Central limit theorems for eigenvalues of deformations of 

Wigner matrices* 

M. Capitaine^, C. Donati-Martin^ and D. Feral^ 

Abstract 

In this paper, we explain the dependance of the fluctuations of the largest eigenvalues of a 
Deformed Wigner model with respect to the eigenvectors of the perturbation matrix. We exhibit 
quite general situations that will give rise to universality or non universality of the fluctuations. 

1 Introduction 

In a previous paper [C-D-Fj . we have studied the a.s. behaviour of extremal eigenvalues of finite rank 
deformation of Wigner matrices and in the particular case of a rank one diagonal deformation whose 
non-null eigenvalue is large enough, we established a central limit theorem for the largest eigenvalue. 
We exhibit a striking non-universality phenomenon at the fluctuations level. Indeed, we prove that 
the fluctuations of the largest eigenvalue vary with the particular distribution of the entries of the 
Wigner matrix. Let us recall these results. The random matrices under study are complex Hermitian 
(or real symmetric) matrices (Mjv)Ar defined on a probability space (SI, J^, P) such that 

Mjv = ^+A^ (1.1) 



where the matrices Wjv and Ajv are defined as follows: 

(i) Wjv is a X Wigner Hermitian (resp. symmetric) matrix such that the N'^ random variables 
{Wn)u, \/23fie((VFjv)jj)i<j, V2'^m{{WN)ij)i<j (resp. the ^'■^'''^^ random variables -^{WN)ii, 
(HOv)ij, i < j) are independent identically distributed with a symmetric distribution fi of vari- 
ance fj^ and satisfying a Poincare inequality; 

(ii) Aat is a determiiristic Hermitian (resp. symmetric) matrix of fixed finite rank r and built from 
a family of J fixed real numbers 6i > ■ ■ ■ > 9j independent of N with some jo such that 
9jg = 0. We assume that the non- null eigenvalues 9j oi An are of fixed multiplicity kj (with 
'l2j^jo ~ ^) A AT is similar to the diagonal matrix 

diag(6'i/fci, . . . , 0jo-i^fejo-i' ^Jo+i^fejo+i' ■ ■ ■ ^^J^kj)- (1-2) 



The Poincare inequality assumption was fundamental in the approach of |C-D-F| . In the present paper, 
we only rely on the results of [C-D-Fj without making use of the Poincare inequality. Hence, we refer 
the reader to jC-D-F] and the references therein for details on such an inequality. Nevertheless, note 
that the Poincare inequality assumption implies that fi has moments of any order (cf. Corollary 3.2 
and Proposition 1.10 in |L]) and this last property will be used later on. 

In the following, given an arbitrary Hermitian matrix M of size N, we will denote by Ai(M) > • • • > 
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Xn{M) its N ordered eigenvalues. 

As the rank of the Ajv's is assumed to be finite, the Wigner Theorem is stih satisfied for the Deformed 
Wigner model (Mjv)Ar (cf. Lemma 2.2 of 0). Thus, as in the classical Wigner model (A^r = 0), the 
spectral measure J2iLi '^Ai(Mjv) of M^r converges a.s. towards the semicircle law whose density 
is given by 

^(x) = -l^V4^^^1[_2.,2.](x). (1.3) 

Nevertheless, the asymptotic behavior of the extremal eigenvalues may be affected by the perturbation 
Ajy. When Ajy = 0, it is well-known that the first largest (resp. last smallest) eigenvalues of the 
rescaled Wigner matrix Wjv/-\/iV tend almost surely to the right(resp. left)-endpoint 2a (resp. —2a) 
of the semicircle support (cf. [B]). This result fails when some of the 0-,'s are sufficiently far from 
zero. Define 

2 

Po, = 0, + ^. (1.4) 

Observe that pg. > 2a (resp. < —2a) when 9j > a (resp. < — cr) (and pg. = ±2a if 9j = ±a). 

For definiteness, we set ki + ■ ■ ■ + kj_i := ii j ~ 1. In [C-D-Fj . we have established the following 

universal convergence result. 

Theorem 1.1. (A.s. behaviour^ Let J+a (resp. J-a) be the number of j's such that 9j > a (resp. 

Oj < -<y)- 

(1) VI < j < J+0-, VI < i < fcj, Afci+...+fc^_i+i(Mjv) — > pg^ a.s. 

(2) Afc,+...+fc,^^+i(Mjv) — >2a a.s. 

(3) \k,+-+kj^j_^iMN) — >-2a a.s. 

(4) Vj > J - J-a + 1, VI < i < kj, \k^+...+k,-i+i(M.N) — > P0J a.s. 

In the particular case of the rank one diagonal deformation Ajv = diag(6', 0, • • • ,0) such that 6 > a, 
we investigated the fiuctuations of the largest eigenvalue of Mjv (with Wat satisfying (i)) around its 
limit pg. We obtained the following non- universality result. 

Theorem 1.2. (CUT) Let Ajv — diag(6',0, • • • ,0) and assume that 6 > a. Define 

^ t /m4-3a\ t a^ 

where t — A (resp. t — 2) when W^r is real (resp. complex) and :— J x^dfj,{x). Then 

cgVN(^Xi{MM)- pg) ^ {M*AA(0,we)}. (1.6) 

Remark 1.1. The strong assumption on the distribution fj, (Poincare inequality) of the entries of 
Wtv is a technical assumption we needed to prove the a.s. result, Theorem M.lY we conjecture it is true 
under more general assumptions, cf. \C-D-F^ ) but the proof of the fluctuations of Theorem \1.2\ only 
requires standard assumptions (existence of the fourth moment) once we know the a.s. convergence. 

On the other hand, in collaboration with S. Peche, the third author of the present article has stated in 
[Fe-Pe| the universality of the fluctuations of some Deformed Wigner models under a full deformation 
A^v defined by (Ajv)ij — d/N for all 1 < i, j < N (see also [Fu-K ). Thus in the non-Gaussian setting, 
the fluctuations of the largest eigenvalue depend, not only on the spectrum of the deformation A^v, 
but also on the particular definition of the matrix Ajv- 

In this paper, we try to explain this dependance of the fluctuations of the largest eigenvalues of 
the Deformed Wigner model Mjv with respect to the eigenvectors of the matrix Ajv- We investigate 
two quite general situations for which we exhibit a phenomenon of different nature. 
First, when the eigenvectors associated to one of the largest eigenvalues of Ajv, say 9j > a, are 
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not "spread" namely belong to a subspace generated by a fixed number Kj (independent of N) of 
canonical vectors of and are independent of iV, we establish that the limiting distribution in the 

fluctuations of Afc^-| |_fc^-_j_|_i(]V[jv), 1 < i < A:^, around pg. is not universal and we give it explicitely 

in terms of these eigenvectors and of the distribution of the entries of the Wigner matrix. 

Secondly, if Kj defined above depends on iV, if there is no "leading" direction among the eigenvectors 

associated to 9j, we establish the universality of the fluctuations of A^^^ hk -i+i{^N), 1 < « < kj. 

We detail these results in the following section. Actually, we assume that the eigenvectors associated 
to the largest eigenvalues 6*1, ... , Oj^^ of Ajv belong to a subspace generated by k{— k{N)) canonical 
vectors of C^. In our approach, we need to isolate aN— kxN — k Deformed Wigner matrix M^-k 
where the eigenvalues of the perturbation are all smaller than cr; we use several well known limiting 
results when N — k tends to infinity involving AlN-k- Hence, our study does not include the full 
deformation case of |Fe-Pe| where k = N. Moreover for technical reasons we have to assume that 
k <C but we conjecture that our result still holds if k N. 

The same kind of questions has been previously studied for the spiked population models by 
[B- Ya2) . The Deformed Wigner matrix model may be seen as the additive analogue of the spiked 
population models. These are random sample covariance matrices {Sn)n defined by Sn — jfY^Yjy 
where Yat is a p x N complex (resp. real) matrix (with N and p = pN of the same order as A'^ — > oo) 
whose entries satisfy first four moments conditions; the sample column vectors are assumed to be i.i.d, 
centered and of covariance matrix a deterministic Hermitian (resp. symmetric) matrix Ep having 
all but finitely many eigenvalues equal to one. The analogue of Theorem 11.11 was established by J. 
Baik and J. Silverstein in [Bk-Slj : when some eigenvalues of Ep are far from one, the corresponding 
largest eigenvalues of Sn a.s. split from the limiting Marchenko-Pastur support. Fluctuations of 
the eigenvalues that jump have been recently found by Z. Bai and J. F. Yao in |B- Ya2| : the setting 
considered in [B-Ya2 ] is the multiplicative analogue of the particular case "fc finite independent of Af" 
in our Theorem 12.11 note that they exhibit universal fluctuations (we refer the reader to [B-Ya2| for 
the precise restrictions made on the definition of the covariance matrix Ep). 

Note that the first steps of our approach as well as the approach of [B-Ya2j are in a spirit close to the 
one of [P] and |B-B-P| . 

The paper is organized as follows. In Section 2, we present the main results of this paper and give 
a summary of our approach. In Section 3, we introduce preliminary lemmas and fundamental results 
which will be of basic use later on. Section 4 is devoted to the proof of Theorem 12.11 and Theorem 
12.21 Finally, we prove some technical results in an Appendix. All along the paper, the parameter t is 
such that t — A (resp. < = 2) in the real (resp. complex) setting and we let :— J x^dp.{x). 

2 Main results 

As in Theorem ll.il we denote by J+o- (resp. J^a) the number of j's such that dj > a (resp. 9j < —a). 
Set fc+CT := fci + • • • + kj^^. We also denote by (e^; i = 1, . . . , A^) the canonical basis of C^. 

We introduce k > fc+o- as the minimal number of canonical vectors of needed to express all the 
eigenvectors associated to the largest eigenvalues 9i, . . . , Oj^^ of Ajv- Without loss of generality, we 
can assume that these fc+o- eigenvectors belong to Vect(ei, . . . , Cfc). This follows from the invariance 
of the distribution of the Wigner matrix W^r by conjugation by a permutation matrix. 

All along the paper we assume that k <C VN- Let us now fix j such that 1 < j < J+a ■ We shall 
study two cases: 

Case a) the orthonormal eigenvectors vf, 1 < i < kj, oi An associated to dj depend on a finite 
number Kj (independent of N) of canonical vectors among (ei, . . . , Cfc) and their coordinates are in- 
dependent of A^ ("The eigenvectors don't spread out"). Without loss of generality, we can assume 
that the vf , I < i < kj, belong to Vect(ei, . . . , exj)', 

Case b) the orthonormal eigenvectors vf, 1 < i < kj, belong to Vect(ei, . . . , e^f^) where Kj = 
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Kj{N) oo when N ^ oo and the coordinates satisfy: 

Vi < kj,\fl < Kj, \{vl,ei)\ ^ as iV ^ oo. 

( "There is no leading direction among the eigenvectors" ) . 

Therefore, we assume that there exists a unitary matrix Uk of size k such that 

diag(C/fc,/Ar-fc)AArdiag(C/fc,/Ar_fc) = diag(6'j4^. , (^j/fe, );< ZN^k+„) (2.1) 

where Zjv-fc^^ is an Hermitian matrix with eigenvalues strictly smaller than 0j_^^. 
In the Case b), Uk satisfies 

k K 

maxmax |(?7fc)ip| — >0 as N ^ oo. (2-2) 

p—l i—l 

Considering now the vectors as vectors in C^J , we define the Kj x kj matrix 

UK,.k,:={vi,...,vQ (2.3) 
namely Uk^ x kj is the upper left corner of Uk of size Kj x kj . It satisfies 

U*K,.k,UK,.k,=Ik,. (2.4) 

Example: 

An = dia.g{Ap{9i),92lk2iON-p~k2) 

where Ap{9i) is a matrix of size p defined by Ap{9i)ij — 9i/p, with 9i,92 > a, p Vn. Then 
= p + fe, ki — 1, K\ — p, — k2. For = 1, we are in Case a) if p does not depend of N and in 
Case b) \i p — p{N) — > +oo. For j = 2, we are in Case a). 

From Theorem 1 1.11 for all 1 < z < kj, Xki-\ i-fcj-i+j(MAr) converges to pe^ a.s.. The main results 

of our paper are the following two theorems. Let ce- be defined by 

%=^- (2.5) 

In Case a) (which includes the particular setting of Theorem ll.2p . the fluctuations of the corre- 
sponding correctly rescaled largest eigenvalues of Mat are not universal. 

Theorem 2.1. In Case a): the kj -dimensional vector 

(^ce^ViV(Afei+...+fc^._i+j(MAr) - pg.); i = l,...,kj'j 

converges in distribution to {Xi{VkjXkj)',i = 1,. . .kj) where Xi{VkjXkj) are the ordered eigenvalues of 
the matrix Vk^xkj of size kj defined in the following way. Let be a Wigner matrix of size Kj 

with distribution given by p (of {i)) and Hk^ be a centered Hermitian Gaussian matrix of size Kj 
independent ofWxj with independent entries Hpi, p < I with variance 



, 2\ t/m4~3a'^\ t 
upp - h(Hpp) _ -(^ ^2 ) + 2 6>2 _ ^2 ' P - - '^i' 



vpi = Ei\Hpi\^) - ^rr^> l<P<l<K, 



(2.6) 



Then, Vk^xkj is the kj x kj matrix defined by 

Vk,xk, = f/^,xfc,(W^if. + HK,)UK,xk,. (2.7) 
Case b) exhibits universal fluctuations. 
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Theorem 2.2. In Case b): the kj -dimensional vector 

(c0^VN{\k,+...+k,-i+i{MN) - p0-);i = l,.. .,kj^ 

converges in distribution to (XiiVkjXkj)', i = l,...fcj) where the matrix Vk^xkj is distributed as the 
GU(0)E(k,xk„fi^). 

Remark 2.1. The condition 1 <C fc <C VN is just a technical condition and we conjecture that our 
result still holds i/ 1 <C fc ^ -/V. 

Remark 2.2. Note that since fi is symmetric, analogue results can be deduced from Theorem \2.1\ and 
Theorem \2.B dealing with the lowest eigenvalues o/Mjv and the 6j such that 9j < —a . 

Before we proceed to the proof of Theorems 12.11 and 12.21 let us give the sketch of our approach 
which are similar in both cases. To this aim, we define for any random variable A, 

ejv(A) =C9^VA^(A-pe,) (2.8) 

with Co- given by (j2.5p . We also set kj-\ := fci + . . . + fcj_i with the convention that fco = 0. 
The reasoning made in the setting of Theorem 11.21 (for which k — k^^ = 1) relies (following ideas 
previously developed in [P] and [B-B-P| ) on the writing of the rescaled eigenvalue ^n{M(J^n)) in 
terms of the resolvent of an underlying non-Deformed Wigner matrix. The conclusion then essentially 
follows from a CLT on random sesquilinear forms established by J. Baik and J. Silverstein in the 
Appendix of [C-D-F| (which corresponds to the following Theorem 13.21 in the scalar case). In the 
general case, to prove the convergence in distribution of the vector {(.Ni^i.. j+i(MAf));* = li ■ • ■ 7%)i 
we will extend, as [B- Ya2] . the previous approach in the following sense. We will show that each of 
these rescaled eigenvalues is an eigenvalue of a kj x kj random matrix which may be expressed in 
terms of the resolvent ofaN — kxN — k Deformed Wigner matrix whose eigenvalues do not jump 
asymptotically outside [—2(7; 2a]; then, the matrix Vk^xkj will arise from a multidimensional CLT on 
random sesquilinear forms. Nevertheless, due to the multidimensional situation to be considered now, 
additional considerations are required. Let us give more details. 

Consider an arbitrary random variable A which converges in probability towards pe ■ Then, applying 
factorizations of type ()3.ip . we prove that A is an eigenvalue of Mjv iff (.nW is (on some event having 
probability going to 1 as ^ oo) an eigenvalue of a kj x kj matrix Afc^^jv(A) of the form 

Xk,,N{X) - Vk^M + i?fc,,iv(A) (2.9) 

where Vk^^N converges in distribution towards Vkjxkj and the remaining term Rkj,NiX) turns out to 
be negligible. Now, when kj > 1, since the matrix Xkj,N{X) (in (|2.9p ) depends on A, the previous 
reasoning with A = A^, ^^.(Mat) for any 1 < i < kj does not allow us to readily deduce that the kj 
normalized eigenvalues ^7v(A^. ^^-(Mat)), 1 < i < kj are eigenvalues of a same matrix of the form 
Vkj,N + op(l) and then that 

(eA'(Afc^_^+,(Mjv)); 1 < * < fc,) = i^^iVk,,N)] 1 < * < fcj) + op(l). (2.10) 

Note that the authors do not develop this difficulty in (B-Ya2| (pp. 464-465). Hence, in the last step 
of the proof (Step 4 in Section 4), we detail the additional arguments which are needed to get (|2.10|) 
when kj > 1. 

Our approach will cover Cases a) and b) and we will handle both cases once this will be possible. 
In fact, the main difference appears in the proof of the convergence in distribution of the matrix V^^^jv 
which gives rise to the " occurrence or non-occurrence" of the distribution fi in the limiting fluctuations 
and then justifies the non-universality (resp. universality) in Case a) (resp. b)). 

The proof is organized in four steps as follows. In Steps 1 and 2, we explain how to obtain 
p.9P : we exhibit the matrix Xkj,N and bring its leading term Vk^^N to light in Step 2. We establish 
the convergence in distribution of the matrix Vk^.N in Step 3. Step 4 is devoted to the concluding 
arguments of the proof. 
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3 Basic tools 



In this section, we fix some notations and recall some basic facts on matrices and some results on 
random sesquilinear forms needed for the proofs of Theorems 12.11 and 12.21 

3.1 Linear algebra 

For any matrix M G Ai]s^{C), we denote by Tr (resp. tr^r) the classical (resp. normalized) trace. 
||M|| is the operator norm of M and ||M||//5 := (Tr(MM*))i/2 the Hilbert-Schmidt norm. Spect(M) 
denotes the spectrum of M. 

For z £ C\Spect(M), we denote by Gm{z) = {zIn — M)^^ the resolvent of M (we suppress the index 
M when there is no confusion). 

Lemma 3.1. Let M be an Hermitian matrix and x G M such that x > Xi{M); we have 

\\G{x)\\ < ^ 



x~ \i{M) 

For Hermitian matrices, denoting by the decreasing ordered eigenvalues, we have the Weyl's 
inequalities: 

Lemma 3.2. (cf. Theorem 4-3.7 of \H-J^ ) Let B and C he two N x N Hermitian matrices. For any 
pair of integers j, k such that 1 < j, k < N and j + k < N + 1, we have 

Xj+k-i{B + C) < X,{B) + Xk{C). 

For any pair of integers j, k such that 1 < j, k < N and j + k > N + 1, we have 

\j{B) + \k{C) < \j+k-N{B + C). 

In the computation of determinants, we shall use the following formula. 

Lemma 3.3. (cf. Theorem 11.3 page 330 in IB-S0j ) Let A G Alfc(C) and D be a nonsingular matrix 
of order N — k. Let also B and *C be two matrices of size k x (N — k). Then 

det ( ^ n ) = dct(i:') dct(A - BD-^C). (3.1) 



D ^ 

3.2 Results on random sesquilinear forms 

In the following, a complex random variable x will be said standardized if E(x) = and E(|xp) = 1. 

Theorem 3.1. (Lemma 2.7 \B-S1^ ) Let B = (6.^) be a Nx N Hermitian matrix and Yn be a vector of 
size N which contains i.i.d standardized entries with bounded fourth moment. Then there is a constant 
K > such that 

E\Y;jBYn - TrSp < KTt{BB*). 

This theorem is still valid if the i.i.d standardized coordinates Y(i) of Y/v have a distribution depending 
on N such that supjv E(|y(i)|'^) < oo. 

Theorem 3.2. (cf. W- YaSS^ or Appendix by J. Baik and J. Silverstein in JC-D-F^ in the scalar case) 
Let A — (ttij) be a N X N Hermitian matrix and {{xi, yi),i < N} a sequence of i.i.d centered vectors 
in X C^^ with finite fourth moment. We write Xi = (xu) G C"'^ and X{1) = {xn, . . . ,xin)'^ for 
1 < I < K and a similar definition for the vectors {Y{1), I < I < K}. Set p{l) — K[xiiyii]. Assume 
that the following limits exist: 



(i) uj = limAT^oo ^ Ej=i o-h 



N I 12 



(a) 9 = limAT^oo :^TrA2 = liniAr^oo jj 1"-^ i ^ 

(Hi) T = liniAT^oc- jjTi-AA'^ = liniAr^oo Y^f,j=i '^iy 
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Then the K -dimensional random vector (l)* AY (l) — p{l)TiAj converges in distribution to a 

Gaussian complex-valued vector G with mean zero. The Laplace transform of G is given by 

Vc e C-^, E[exp(c^G)] = cxp(ic^Bc), 

where the K x K matrix B = {B{1, V)) is given by B = Bi + B2 + B3 with: 

Bi{l,l') = uj{E[xnynxinyi'i]- p{l)p{l')) 

B2{l,l') = (e-Lu)E[xnyi,i]E[xinyii] (3.2) 
B3{l,l') = (t - Lj)E[xnXin]E[ynyi'i]- 



4 Proofs of Theorem 12.11 and Theorem 12.21 

As far as possible, we handle both the proofs of Theorem 12.11 and Theorem 12.21 We will proceed in 
four steps. First, let us introduce a few notations. 

For a TO X g matrix B (or B) and some integers 1 < p < m and 1 < ^ < 97 we denote respectively by 
[-^]pxJ' i^]pxi^ i^]p'xi ^^'^ [^]pxi upper left, upper right, lower left and lower right corner of size 
p X / of the matrix B. If p = /, we will often replace the indices p x I by p for convenience. Moreover 
if p = TO , wc may replace or \ by — )■ and or \ by <— . Similarly if ^ = g, we may replace or 
\ by t and ^ or \ by i. 

For simplicity in the writing we will define the k x k, resp. N — k x N — k, resp. k x N — k matrix 
Wk, resp. WN-k, resp. Y, by setting 

w^^fllSl (41) 



, Y* WN-k , 

Given B e Aljv(C), we will denote by B the N x N matrix given by 

B diag([/,*,/^_fc)Bdiag(t/,,/^_fc) ^( f l'^^^^"^' 

V ^N-kxk iJN-k 

One obviously has that B^-k = B^-k- 

In this way, we define the matrices Mjv, Wjv and Ajy. In particular, we notice from (|2.ip that 



An = dmg{e,Ik^,{9ak,)i<j,.,i^,,ZN-k^^) = ( f 'I'-''-' ). (4.2) 

V ^N-kxk AN-k J 

Note also that since Ajq-k is a submatrix of ZN-k+^j all its eigenvalues are strictly smaller than a. 
Let < 6 < [pq. — 2f7)/2. For any random variable A, define the events 

^%\\) = |ai (^^+diag(C/fe,/Ar_,)diag(Ofe^„,ZAr_fc^Jdiag(C/,*,/Ar_fc)^ <2a + 5-\> pg^ - 5^. 
= [x,{^+AN-k)<2a-V6 



and 



f]w(A) = r!j^'(A)f|f)j^'. (4.3) 



On fiAr(A), neither A nor pg. are eigenvalues of M^-k '■— + ^N-k and the resolvents G(A) := 

{Xl^-k — M^-k)^^ and G{pgj) :— {pg^lN-k — M^-k)^^ are well defined. Note that from Theorem 
11.11 for any random sequence A^r converging towards pg- in probability, limjv >oo IP(f^Af(AAr)) = 1. 
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STEP 1: Let A be a random variable. On ilN(X), 
det(Mjv - A/at) = det(MAr - A/at) 



det 



det{MN-k - A/jv-fe) det (m^ - Xh + AIkxN-kG{X)MN_kxk) ■ 



The last equality in the above equation follows from p.ip . Since on Qn{X), X is not an eigenvalue of 
M^^k, we can deduce that A is an eigenvalue of Mjv if and only if it is an eigenvalue of 



Qfc,Ar(A) := Mk + MkxN-kG{X)MN-kxk- 

Now, note that we have also from (|3.ip that 

+ ^Ar-fe+„ — XlN-k^ 



(4.4) 



det 



W 



N 



det 



N 

WN-k 



N-~k+, 



N 



[ZN-k+^]^_k - A/AT-fc^ X det ([Qfe.7v(A)] 



Xlk- 



The matrix 



N-k + a 



ZN-k+„ is a submatrix of +diag(Ofe^^, ZN^k+„) whose eigenvalues are 



(on r2jv(A)) smaller than 2a + S. So, since on CIn{X), X is greater than pg, ~5 > 2(7 + 5, we can conclude 

... N 

that A cannot be an eigenvalue of 
we can define 



^ ^ +Zjv_fc+„, and then neither of [Qfe^jv(A)]^^.^^. Thus, 



Efc_fc^^(A) := ([Qfc,Ar(A)]>l, -A4 



(4.5) 



Moreover on f7jv(A), one can see using p.ip that if Aq is an eigenvalue of [Qfe,Af (A)]^/!^^ — XIk-k+„ 
then A is an eigenvalue of 



W 



N 



N-k+, 



Hence, 

and then 

so that finally 



+ ZN_k+^ - diag(Ao/fe-fc+„,07v-fe)- 



Ao \ > POj - 5 - 2a - 6, 
1 



.(A)|| < 



(4.6) 



pg^ -2a -26 

Using oncemore (I3.ip . we get that on r2jv(A), A is an eigenvalue of Qfe.Ar(A) if and only if it is 
an eigenvalue of [Qfc,7v(A)]^^ - [Qfc,Af (A)]^^xfe-fe+„ ^fc~fc+„(A) [Qk,NW]'{_k^^-^k+„ or equivalently if 
and only if ^Ar(A) is an eigenvalue of 

CO,VN([QkM^)]>-PO,h,. ~ [QkMW]Cxk-k^. ^k-k^X) [Qk.N{X)]C 



Now using 

G(A) - G{pe^) = -{X- peMpe,)G{X), 

one can replace G(A) by G(peJ + -(A - pej)Gipe^) {g{pbj) - (A - pe.)G{pe.)GiX) 
following writing 

1 1 „ N-k cr2 



and get the 



^YGiX)Y*^^YGipe^)Y*+^.iX)DkM^)-^.iX)- ^ ^ 



(4.7) 
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where 



Then 



1 



92 - (t2 



VI (Wu + ^{YG{pe,)Y* -{N- k)^h 
+VNdiag (Ofe^ AOi-0j)hnl = h---.J+.J^ j) 



eiv(A). 



N 



cr^ ce, N 



92 - a2 



The following proposition (adding an extra matrix A^^^ for future computations) readily follows: 
Proposition 4.1. For any random variable A and any fc+o- x k+a random matrix Afc^^, on CIn{X), 



A is an eigenvalue o/M^r + diag(Afc^^ , 0) i/f C-/v(A) is an eigenvalue ofli-k 



.N 



(A) 



NAk^„ where 



Xfc^„,A,(A) := [t/;Sfc,^t/fe]^ + ViVdiag (Ofc^., (0; - 9,)h, , i ^ I, . . . , J+^J ^ j) + ^nW [U*DkMWUk_ 



+ 



{ef- 



where 



B. 



k,N 



Wk 



-^f-j Ik+„ — :^rfe^^xfc-fc+<,(A)Sfe-fe_|_„(A)rfe^^xfe-fe+„(A)* 

(4.9) 
(4.10) 



(4.8) 



^YG{pe^)Y*-{N-k)-h 



ce^Dk^NiX) = Tw(A) + 0Ar + V'jv 
TNiX) ^ ^{X^ Pe,)YGiX)G{pe^rY* 



9 AT 



rG(peJ'l'*-^'TrG(peJ'4 



N 



2N-kf^ ^. ,2 1 



N 



tTN-kG{p0^) - -J 



9 y ^fei 



and 



^k+^xk-k+AX) = Tn{X) + Ak^^{X) with 



(4.11) 



Tat (A) 



Ul{Wk + -^YG{X)Y*)Uk 



Afc+„(A) = [/,!yG(A)iAr_fcxfc 



Moreover, the k — fc+o- x fc — fc+o- matrix "Ek-k+^iX) defined by (|4.5p is swc/i f/iai 

||Sfc-fc+„(A)|| <l/(pe, -2a-2(5). 

Let us make some comments on our approach in order to explain why we proceed in two steps 
namely we apply twice a factorization of type p.ip to deal with a fc x fc matrix and then with a 
fc+o- X fc+o- matrix. This approach makes the accommodating resolvent of the Deformed Wigner matrix 
+AN-k arise. A single application of a factorization of type (|3.ip to go from a x iV to k+a x k+a 
matrix would require to deal with the matrix [^N/VN]^_f,^ + Z^-k^^ whose limiting spectral 

and 



behaviour is a priori unknown. Moreover, the independance of the matrix [Wn/\/N]^ 



cx JV-fc 
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— j^+Api-k arising from our two steps approach will be of fundamental use in the following, whereas 

their analogues in a single step approach are [WAr/-\/]V]fc^^xW-fc+„ ^^'^ ['^n/VN]^_i^_^^ + Z^-k^^ 
which are not independant. 

Throughout Steps 2 and 3, Ajv denotes any random sequence converging in probability towards 
pe-. The aim of these two steps is to study the limiting behavior of the matrix Xfc^^jv(AAr) as N 
goes to infinity. 

STEP 2: We first focus on the negligible terms in Xfc^^_jv(AAr) and establish the following. 

Proposition 4.2. Assume that k <^ \/N . For any random sequence A^r converging in probability 
towards pg., on f2Ar(Ajv), 

Xfc^„,jv(Aw) = Vk^^.N + %/iVdiag {Ok.A^i ~ Oo)h, , Z = 1, . . . , Z ^ j) + (1 + \S,n {^N)\f or{l) . 

(4.12) 

with Vk^^^N given by 

Vk^^,N ■■= [U^BkM]',:^^ ■ (4.13) 

The proof of this proposition is quite long and is divided in several lemmas. Although our final 
result in the case k infinite holds only for k ^ VN, we will give some estimates for fc ^ iV once this 
is possible. 

Lemma 4.1. Let k ^ N. Then, on ^^{An), 

[U*,Dk,N{^N)Uk]'^_,^ = op(l). (4.14) 

Proof of Lemma 14. It We refer to Proposition 14. II for the definition of Dk,N{AN), tn, <Pn and tpN- 
Let K = diag(/fc+^,Ofc-fe+J. On f7Ar(AAr), 

II K^iv^7fc]^^^ Whs = ^IAat - peMTT{KU:YG{AN)G{pe,)'Y*UkKUlYG{pe,)'G(AN)Y*UkK)}^ 

< l|Ajv-peJI|G(pejin|G(AAr)||||y*c/feif{/,*y||^{Tr(/^t/fc*yr*c/fci^)}^ 

< 1^\An - peMG{peM'\\Gi^N)\\TT{KUlYY*UkK) 

< ^|A--P«.l(^^-^32^Tr(i.f/*yr*L/.X). 

We have, 

^ ^ N-k fc+<7 k 

-TT{KUlYY*UkK) = - ^ ^ 5] (t^W^/,fc+pW^ (C/fe),,. 

k+a ^ N~k k 

= ^J^Y^ H {'Uk)l,iWl,k+pWq^k+p {Uk)q,i 
i—1 p—1 l,q—l,l^q 

N~k k 

+E^EEi(^'^)Mpr..+pP- 

i=i p=i /=i 

Since {X]f g=i i^q {Uk)i,iWi^k+pWq^k+p {Uk)q.i] 1 < p < N — k} are i.i.d random variables with mean 
zero and such that the second moments are bounded in N , we can deduce by the law of large numbers 
that jfY^pZi Y^\^q^ij^q (Uk)i,iWi^k+pWq.k+p {Uk)q,i couverges in to zero and thus in probability. 
Similarly, since {X]f=i l('^fe)i,iP|W'i,fc+pP; 1 < p < N — k} are i.i.d random variables with mean 
and bounded (in N) second moments, by the law of large numbers Y^p^i X]/Li l(C^fc)i.iP|l^i,fc+pP 
converges in towards and thus in probability. It follows that when N — > +cx3, 

^Tr{KUlYY*UkK) ^ k+„a^ (4.15) 
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Hence [UlTNUk]^^^ = op(1). 

It follows from Lemma |5. II in the Appendix that 

,N-kr c,, 1 



[Uk^NUk]^^^ -a -^^^tTN-kG{pe^) - ^j^-^J = op(1). 



Now, we have 

E 



(ll [U*k<l>NUk]'^_^^ ln„(A„)||lf5) < E (ll mc^NUk]^^^ l<lll?s) 

fc+. 1 / 



k , „ 

where for any p = l,...,k+cr, we let U{p) = *[{Y*Uk)i,p, ■ ■ ■ , {Y*Uk)N-k,p]- We first state some 
properties of the vectors U (p) . 

Lemma 4.2. Let hi denote the N — k x k-^-^r matrix [Y*Uk]kr^ ■ Then, the rows {Ui,;i < N — k) 
are centered i.i.d vectors in C'^+'', with a distribution depending on N. Moreover, we have for all 
1 <p,q < k+a- ■■ 



E{UipUiq) = Sp ga"^ with E{UipUiq) = in the complex case, 

E^prnql'] = (1 + yp^q)a^ + [E{\W,,\^) - (1 + i)a*] ^ \{Uk\pnUk)i,q\' . (4.16) 



1=1 



Since X](=i — 1; fourth moment oflAip is uniformly hounded. 

We skip the proof of this lemma which follows from straightforward computations using the indepen- 
dence of the entries of Y and the fact that Uk is unitary. 
Then, according to Theorem 13.11 and using Lemma 13.11 



^E(|Z^(p)*G(p,j2^(p)-a2TrG(pe,)'plo„(A„)) < ^^{tvr, G{pe,)%,^ 

^E(\\G{pe,)\ 



< '-^E{\\G(p,MH 



K 

< 



N (pg^ -2a- <5)4 • 

Besides for p ^ q, using the independence between {U{p),U{q)) and G{pg.), we have: 

N-k 

E(\U{pyG{pe^)^U{q)\H^^2>^ = E[U,p{G%UjgUipW)^U^gl^i2,] 

N-k 

= Yl H^ipl^jql^lp^mq]E[{G%jW)^l^i2,] 

where we denote by G the matrix G{pg-) for simplicity. ^From Lemma |4.2[ for p ^ q, the only terms 
giving a non null expectation in the above equation are those for which: 

1) i = I, j = m and i j. In this case, 

E[U,pUjgU^pUjg] = E[a^pU^p]E[UjgUjg] ^ CT^ 

and 



N-k 

Y E[{G%,Whin^^^]<ETT{G%,2,). 
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2) i — j = k = I. In this case, using (|4.16p . there is a constant C > such that 



Moreover 

N-k 

^E[G|G?,V,]<ETr(G4X 



i=l 

Therefore, 



E 

Hence, 



(^U{prG{pe^fU{q)\H^^.^) < {C + a^) ETr{G{pe^)%,.>). (4.17) 



Thus 



4\ + 



The convergence in probabihty of [U^cjyNUk]^^ towards zero readily follows by Tchebychev inequality. 
Lemma [4. II is established. □ 



For simplicity, wc now write 

S(Aw) = Sfc_fc^„(Ajv) 

Let us define 



Rk.N (Aat) j=—Ik^^ + --^ '^Ar( n) j^^^ — rfc_|_^ xfc-fc+„ (AAr)i;(Ajv)rfe^^ X fe-fe+„ (Aat)* 

V iV — cr eg- 1\ 

(4.18) 



To get Proposition 14. 2[ it remains to prove that if fc <C v A^, 

Rkn{^n) - (1 + \^N{^N)\for{l). (4.19) 
Once k <C "\/iV, we readily have that 

+ i2 12 ~ = (1 + K^(AJV)I) Op(l). 



Hence, (|4.19p will follow if we prove 

Lemma 4.3. Assume that k ^ ^/N. On ^^.^{An), 

-^rfc+„xfc^fc+„(Ajv)S(A^)rfe+„xfc-fe+„(AAr)* = (1 + |^jv(AAr)|)2op(l). (4.20) 



For the proof, we use the following decomposition (recall the notations of Proposition 14. ip : 



rfc+„xfe-fe+„(AAr)i;(AAr)rfe^^xfc-fc+„(AAr) 
= 1 

where (using (|4.7|) ') 



TatET^ + TatSA^^^ (Aw)* + Afe^„ (Ajv)SAfe^„ (Ajv)* + Afe^„ (Aw)ET^. (4.21) 



U*,{Wk + ^YG{AN)Y*)Uk 



- fc+CT X (fc — fc+„) 

and we replaced S](Ajv) by S. We will prove the following lemma on Tjv. 
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Lemma 4.4. // /c < \/N , 



II [U*kDu,N{^N)Uk]i 



-[- £7 X fc-j- ( 



||h5 = ov{Ni). 



II KBfc,Ar(A^)C/fe]/ ||h5 - ov{m). 

and therefore, for k ^ \/N, 

\\Tn\\ =Op{Ni){l + \^N{AN)\). 

Proof of Lemma 14. 4t To prove (|4.22p . we use the decomposition 
As in the proof of Lemma 14. 1[ we have 

4% kkj^u 1 



(4.22) 
(4.23) 



E 



N {pe^ -2a- 5)^ ' 
SO that, for k <^ N and using Tchebychev inequality, we can deduce that 

II Pk<l>NUk]k^^y,k_k^^ ||hs1o„(a„) = op(l). 
Let K = diag(/fe^^,0) and L = dia,g{0, Ik-k+ J- On fliyiAN), 



HS = ^lAjv - pe^\{TT{KU'^YG{AN)G{pe^)^Y*UkLU*kYGipe^fGiAN)Y*UkK)}i 

< ^|AAr-peJ||G(peJ||2||G(Ajv)||||y*;7feL;7,!y||^{Tr(i^;7,!yr*C/feX)}5 

< ^|Ajv-pe,|||G(peJ||2||G(Ajv)||{Tr(LC/,!ry*{/fcL)}5{Tr(if{/,*rr*;7feif)}^ 



{Tr(LC/fc* yr* c/fcL)} 5 {Tr(i^ ;7fc*ry * [/fci^)}^ 



According to (|4.15p . ^^{Tr(i4rt/^yy*C/feiir)}2 converges in probabihty towards y/k+^a. 

/c N~k k 

-Tr{LUlYY*UkL) - ^ E E E WkhlWi,k+pW^ iUk),,r 

N — k k k 

+^E E Y.\^Uk)i,?m,k+,\ 



=1 i=i 



{Ei=fc+^+i I]i.g=i,i5^g (C^fc);,iW^;,fc+p Wq^k+p {Uk)q,i, I < p < N - k} aiG i.i.d variables with mean zero 
such that 



E 



fe k 

E E W^^i Wq ,k+p 



< C{k - k+„f 



for some constant C . Hence 
/ 



E 



^ E E E {Uk)l,zWl^k+p Wq^k+p {Uk)q,, 



p—1 i—kj^^-\-ll,q—l,l^q 



< G 



J 



(fc - k+^f 
N 



so that the first term in the previous sum converges in L"^ towards and thus in probabihty. 
Moreover, since X^iLfe +i F^i — S?=i l(t^fc)i,iP|l^!,fc+pP i-i-d random variables with mean and 



bounded second moments (in N), there exists some constant C such that 
/ 2 

N-k k fc 

^ ^E E rr^EK^^Mnw^. 



N ^ ^ k-k,^ 
p=i i=k+„+i ^ 1=1 



|2 ^2 

,fe+p| —a 



< 



C_ 
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so that jfY.p=iY.2=k+^+i k-k+„ Sf^i \iUk)Li\'^\Wi^k+p\'^ converges in towards and thus in 
probabihty. It follows that j^Tv{LUlYY*UkL) = Op(fc) and then that under the assumption k <C ViV, 



Note that we also have ior k ^ N , 



\U*krNUk]{^^^^_r,^^ Whs = op(Vfc) = op(iV3). 
PlrNUk\{^^^^_^^^ Whs = |^w(A)|op(1), 



and therefore 

II [U*kDk^N{t^N)Uk]C^^^_. Whs = (1 + |eiv(A^)|)op(l). 



Thus, ((i:^ is established. 

For recall that [UtBkM]{^^^u-k,. = [UlWkUk]{^^,,_,^+^[UtYG{pe,)Y*Uk]{^^^,_,^^ 



Since 

one has that 



E (\mWkUk]{^^^^_^J\ls) < k+.ka^ 



'{\mWkUk]] 



Hence, as fc < V^, we can deduce that W[UkWkUk]'k_^_^xk-k+J\HS = op{Ni). 

Now, let us prove the same estimate for the remaining term. Using the same proof as in l4.171 one can 
get that for p ^ q, for some constant C > 0, 



E 



(|Zi(p)*G(pe,)Z^(g)pl^(.)) < CETr(G(p9,)'l^(2)) 



and then that for some constant C > 
1 



Ef 



^mYG{po^)Y*Uk]^^^^,_,^^ < Ckk^^ (^^^ _ 2, _ s)^ ■ 

Then using that 

\^[U:YG{pe^)Y*Uk]C^^,_,J\Hsi^^^^>em'^ < -2_E[W^[UlYG{pe^)Y*Uk]{^^^,_,J\ls] 



we deduce since k <ti \/ N that 



Thus (|4.23p and Lemma are proved. □ 
Using that 



W^[UtYG{pe,)Y*Uk]{^^^^_,J\Hs\^i^^ = or{m). 



< 1 -, (4.24) 

Pe, - lu - lo 



one can readily notice that Lemma 14.41 leads to 

1 



^TjvST^ = (1 + |eAr(A^)|)'op(l). (4.25) 



We now consider the remaining terms in the r.h.s of (|4.2ip . We first show the following result where 
we recall that t^k^^Po,) = \U"UG{pe^)Aj^-kxk\C^^^k-k^„- 

Lemma 4.5. ^T7vi;Afe^^(p(,J*, ^Afc^^(pej)SAfc^^(p0j* and ^ A^^^ (pejET;^ are all equal to 
some (1 + |Cjv(Aiv)|)op(l). 
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Proof of Lemma 14.51 : We will show that, on ^^{An), for any u > 0, 

Afe+.(pe,) =op(iV"). (4.26) 

One can readily see that this leads to the announced result combining Lemma [44l (|4.24p and (|4.26p . 
First, using the fact that U^Y is independent of i^(2>G{pQj) and that for any p, the random vector 

U{p) [{Y*Uk)i,p, ■ • ■ , {Y*Uk)N-k,p] has independent centered entries with variance a^, one has that 
E(lo„(A„)TrAfc^„(pe,)Afc+„(pe,)*) < m^iv Tr Afc^^(p,JAfc^„(pe,)*) 

< k+^a^E \^1^^(2) \\G{pe^)fTTAN-kxk-k+A*N~kxk^k^, } 

< IZ Tr 4^ 



^2 



J 



ipe^ -2a- 5Y f^^ 

Therefore, P(1o„(a„)|| Afc^^(pe,)||ffS > eN^) < e~^N-^"E{ln^^A^j\\Ak^^ipe^)\\ls) goes to zero as 
N tends to infinity. Hence (|4.26p holds true on ^^{An) and the proof of Lemma is complete. □ 

Let us now prove that 
Lemma 4.6. Afc+„(Ajv) = Ak^^{pe^) + Op{\^n{An)\). 
Proof of Lemma 14. 6t We have 

AkUAN)-A,^M = ~{^N~Pe,)[U^YG{pe^)G{AN)AN-kxk]C^^k-k+^. 

Let us define Vfc^^ — [C/^yG'(pej )G'(Aiv)^7V-fexfc]fc^^ xfc-fc+<, • Then for some constant C > depend- 
ing on the matrix ^jv-fcxfc, 



Tr(Vfe^„V^^J < C\\Gipe^)r\\GiAN)rTTiU*U) 

C_ 

2^- 2sy 



where we denote as before U = [5^*t/fc]fe^^- Thus letting C" := C Cg^ , 

II Afc,„ {An) - A,.^„ {p6^ )\\ls<C' {Cn (An))' (^^^ _2a-26r ^ ' 



It follows from KTB that 



implying that ||Afc^^(AAr) - Ak^^{pe.)\\HS = Op(|Cjv(Aw)|)- □ 

We are now in position to conclude the proof of Lemma 14.31 Indeed, writing 

Afc^„ (A^)ET^ = (Afc^„ {An) - A,.^^ {pe^ )) ET^ + Afe^^ (pe, )ST^ 



and 



Afc^^(A^)SAfe^^(AAr) = Ak^SPe,)^^k^APe,r 

+ (Afc^„(A^) - Ak^^ipe^)) ^Ak^Ape,)* 

+ (Afc+„(Aiv) - Afc+.(P9,)) S (Afc^„(A^) - Ak^APe,)y 

+ Afc+„ (pe, )S ( Afc^„ (Aat) - Afc^„ (pe, )) * , 
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we deduce from LemmasimiUland (glMl), ((i:^ that ^ Afc^^ (AArjET;;^ and ^ Afc_^^ (A7v)SAfe^„ (Aat)* 
are both equal to some (1 + |^Af(AAr)|)op(l). Using also (I4.25p . we can deduce that 



^rfc^„xfc-fc+„(Ajv)Srfe^„xfe-fe+„(Ajv)* = (1 + |eiv(A^v)|)'op(l) 
which gives (|4.20p and completes the proof of Lemma 14.31 □ 



(4.27) 



Combining all the preceding, we have established Proposition |4?2l We now prove that provided it 
converges in distribution, with a probability going to one as N goes to infinity, ^Ar(Ajv) is actually an 
eigenvalue of a matrix of size kj . 



Lemma 4.7. For all u > 0, 



\Vk+„,N\\HS 



op(l) 



Proof: Straightforward computations lead to the existence of some constant C such that 

E(\\MWkUk].\\Hs) <C. 



The convergence of || [U'^WkUk]).^ 11/-^" in probability towards zero readily follows by Tchebychev 
inequality. Following the proof in Lemma H?T] of the convergence in probabihty of [t^fe^ArC^fc]^,^ 
towards zero, one can get that 



E 



1 



On(An) 



< 



and the convergence in probability towards zero of the term inside the above expectation follows by 
Tchebychev inequality. Since moreover according to Lemma l5.1[ 



1 



(TrG(p«J-(^-fc)-) =op(l). 



we can deduce that 



iV-"|| 



Op(l). 



The proof of Lemma [4. 71 is complete. □ 



Proposition 4.3. Let A^^ be an arbitrary kj x kj random matrix. If^Ni^N) converges in distribution, 
then, with a probability going to one as N goes to infinity, it is an eigenvalue of Xfc^^^7v(AAr) + 
diag(Afc^. ,0) iff ^n{^n) is an eigenvalue of a matrix Xkj,Ni-^N) + A^^ of size kj, satisfying 

Xk,,N{^N)^Vk^^N + oril) (4.28) 

where Vk^^N is the kj x kj element in the block decomposition ofVk^^^N defined by (j4.13p .' namely 

with Ux xk- and B^.n defined respectively by (j2.3p and 



Proof of Proposition 14. 3t Since £,n{-^n) converges in distribution, we can write the matrix 
Xfc^„,Ar(AAr) given by (|4.12p as 



Xfc^„,7v(Aw) = VAfdiag(Ofc^., (((?, - ej)Ik,Wj) + i?fc+„,Ar(Ajv) 
where Rk+„,N{^N) ■= Vk^^^N + op{l). Let us decompose 'K.k+„,N{-^N) in blocks as 



X/c^^,Ar(Ajv) 



a:, 

Xk 



ki,N 



X, 
Xk 



kj X fe+c — fcj ,N 



16 



We first show that £,n{^n) is not an eigenvalue of Xk^„-k.i,N- Let a — inf/^j \9i — 9j\ > 0. Since, 

Xk+„~kj,N = VNdia.g{{{ei - 9j)Ik,)i^j) + Rk+^-k^,N, 
if /i is an eigenvalue of Xk^^-kj, then 

\fi\/VN>a-\\Rk^^.k„N\\/^- 

Now, using Lemma HTTl 



\\Rk^,-k,M\/^N^op{l). 
Hence ^Ni^N) cannot be an eigenvalue oi Xk_^^^kj,N- Therefore, we can define 

^kj,N = Xkj^N — XkjXk+„-kj,N{Xk+„-kj,N ~ ^N{^N)Ik+„-kj) ^ ^k+^-kj xkj ,N 

— Vkj.N ~ RkjXk+„-kj,N{Xk^^-kj,N — ^N{^N)Ik+^-kj) Rk+^-kjXkj ,N + Op(l). 

To get (|4.28p . it remains to show that 

\\RkjXk+a-kj,N{Xk+^~kj,N — £.N{^N)Ik+„-kj) ^ R-k+^-kj xkj ,n\\ ~ Op{l). 

This follows from the previous computations showing that (for some constant C > 0) 

\\{Xk^^-k,M - ^NiAN)Ik+^-k,)-'\\ < {C + op{1))/Vn, 

combined with the definition of Rk+^.Ni^N) and Lemma 14.71 The statement of the proposition then 
follows from (13.111. □ 



STEP 3: We now examine the convergence of the kj x kj matrix Vk-.N = U^.^^.\Bk^N]K-UK xk- in 
the two cases: Kj independent of N and Kj ^ oo. 

a) Kj and the matrix UKjXkj are independent of N 

Proposition 4.4. The Hermitian (resp. symmetric) matrix [Bk^N]^. converges in distribution to- 
wards the law of Wkj + where is a Wigner matrix of size Kj with distribution given by fi 
(cf (i)) and Hk^ is a centered Gaussian Hermitian (resp. symmetric) matrix of size Kj independent 
ofWKj, with independent entries Hpi, p < I with variance 



■ —J " (4.29) 



vpi = n\HpiV) - t;^— ^, l<p<l<K, 



The proof follows from Theorem 13.21 and is omitted. We shall detail the proof of a similar result in 
the infinite case (cf. below the proof of Lemma [49l)- 



b) Kj{^ Kj{N)) ^ 00 and UK,xk,{^ UK,xk,{N)) satisfies O 

Proposition 4.5. Ifms3^pL^ max^^-^ |(J7fe)ip| converges to zero when N goes to infinity then the kj x kj 
matrix U'^ ,^^\Bk^N]x -^KjXkj converges in distribution to a GU(0)E(kj x kj, gi_^2 )■ 

We decompose the proof into the two following lemmas. 

Lemma 4.8. // maXp!^-^ max^^^-^ |(?7fe)ip| converges to zero when N goes to infinity then the kj x kj 
matrix U'^ ,^i^,[Wk]^ UKjXkj converges in distribution to a GU(0)E(kj x kj,a'^). 

Proof of Lemma 14. 8t First we consider the complex case. Let Upq ^C, l<p<q<kj and 
app G M, 1 < p < /cj, and define 

L^ia) ^ iapq{UkWkUk)pq + apq{U^WkUk)pq) + ^ '2.app{UlWkUk)pp. 
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We have 



where 



LN{a) = Y,Di{WN)n+ J2 Rii{V2^e{{WN)u)) + J2 Iii{^'^rn{{WN)ii)), 

i=l l<i<l<Kj l<i<l<Kj 

^<p<q<kj 



l<P<9<fe3 

^ 1 ^^^y^^L-. / 



K ^ 

Hence Lm{ol) = 'Ylm^i l3m,N4'm where (pm are i.i.d random variables with distribution /i and /?m,jv 
are real constants (depending on the Upq) which satisfy max^^^^^j^ l/^m.Tvl ^ when N goes to infinity. 

Therefore the cumulants of Lpfia) are given by C4^^ = J2m=i l^m N^nil^') for any n G N* where 
Cnifj) denotes the n-th cumulant of ^ (all are finite since fj, has moments of any order). In particular 
C^^^ = 0. We are going to prove that the variance of Li^{a) is actually constant, given by 



(N) Kj 

^ = Ea^+ E 4+ E 4 = 2 E E (4-30) 

i=l l<i<j<K-i l<i<j<K-i '^<P<q<kj i<P<kj 

Let us rewrite the l.h.s as 

E^^+ E E ^ = E 



j=l l<i<l<Ki l<i<l<Ki 



P,<I,P ,1 

1 < p < q < kj 



where 

= 2^JRe(a,,([/,),,(C/,),^)5ie(ap.,.(C7fcV(C/,),^.) 

i=l 

+ ^ (ap,((t^fe)i,(C4),p + (UkhMdip)) (t.,-,' ((t/fc),,' (C4),y + (C/fe)^ (C4),y )) 
+ ^ (a^,((C/fc);,(f4)^ - {Ukh.JlWp)) (a^.,. ((«7fe);,' (f^V - (C/fc),,'!^) 

l<i<i<ifj 

So that 

n^,,' = E {3«e(aj,,(([/fe),g(C4),^ + (t/,)i,(t4),^))5Re(ap,,,(([/0;^^ 
i<i,;<Kj 

+9m (ap,((f/fe),,(f4)^ - (f/fc),,(f4)^)) 9m (q^,^, ((C/^c),,' 04)" - (C/fe),,' (t^V )) } 



P,<I,P ,Q P,Q,P ,9 



where 



K1p\,' = E ['^^{^P'i{Uk)i,{Uk\^)^e(a^.^.(Uk\-{Uk)i 

l<i,l<Kj 



l<i,l<Ki 



\'^P<l\'^^{p,q),{p',q') 
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and 

l<iA<Kj 

-3m {apq{Uk)igiUk)rp^ (op',' (?/fc),^' (J/fe)ip' ) } 

^ l<iA<Kj ^ 



Then (|4.30p readily follows. In the following, we let const = J2m=i Pm n- 

Since |C4^^| < consi max^^^^ |/3m,Af|"~^|C'„(/Lt)|, Cn^^ converges to zero for each n > 3. Thus we 
can deduce from Janson's theorem [J] that Lm{q) converges to a centered gaussian distribution with 



variance cr^(2^^^ |apqp + 4^^< |appP) and the proof of Lemma is complete in the 



complex case. 

Dealing with symmetric matrices, one needs to consider the random variable 

LN{a) := ^ a.pq{UiWkUk)pq + ^ app{UlWkUk)pp 



for any real numbers apq, p < q. One can similarly prove that Ln{q) converges to a centered gaussian 
distribution with variance (T^{'iY.i<p<q<k, "pg + 2Ei<p<fc, "pp)- ° 

Remark 4.1. Note that Lemma \4.8\ is true under the assumption of the existence of a fourth moment. 
This can be shown by using a Taylor development of the Fourier transform of (a) . 

Lemma 4.9. // max^^j^ max^^-^ |([/fc)ip| converges to zero when N goes to infinity then the kj x 



kj matrix -^C/^^^^^. 



I. . 2 



(yG{p0.)Y* - {N - k)^Ik^ ^ UKjXk, converges towards a GU(0)E(kj x 



Proof of Lemma 14. 9t We shall apply a slightly modified version of Theorem 13.21 (see Theorem 7.1 
in |B- Ya2| ) but requiring the finiteness of sixth (instead of fourth) moments. Let K — kj{kj + l)/2. 
The set {1, . . . K} is indexed by I = {p, q) with I < p < q < kj, taking the lexicographic order. We 
define a sequence of i.i.d centered vectors {xi,yi)i<N-k in x by xu ~ Uip and yu — Uiq for 
I = {p,q) where U is defined in Lemma 14.21 The matrix A of size — fc is the matrix G{pg ) and is 
independent oiU. Note that we are not exactly in the context of Theorem 7.1 of |B-Ya2| since the i.i.d 
vectors {xi,yi)i depend on N (and should be rather denoted by {xi^M ,yi,N)i) but their distribution 
satisfies: 

1. pit) = WiXiiyii] — Sp^qcr"^ for I — {p,q) is independent of iV. 

2. E[xnyin] = Sp^qXJ^ if I - (p, q)J' - (/, q') (see B2 in 

3. Complex case: E[a;iia;i'i] — ^[ynyi'i] = if / = {p,q)J' — {p',q') (see B3 in (13.21) 1. 
Real case: E[xaXi'i] = cr^f^p^p' and E[y;iy;/i] cr^Sg^g, if I = {p,q),l' = {p' ,q')- 

4. (see Bi in 

' ^[xnynxviyvi] = cr'^{Sp^gSp,^g> + 6p^g,5p,^g)+ 



m\Wi2\^) - 2a^] Y.7=i{Uk\,q{Uk\,p{Uk),^g' {Uk),^p' in the complex case, 

^[xiiyilXl'iyi'l] = a'^iSp^gSp'^q' + Sp^g'Sp'^g + Sp^p'Sq^q') + 



[£(11^121^) - StT^] j:fMUkkq{Ukkp{Uk\q' iUk),^p' in the real case. 
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Under the assumption that max^L^ max^Jj^ \ {Uk)i,p \ converges to zero when k goes to infinity, the last 
term in the r.h.s of the two above equations tends to 0. 

It can be seen that the proof of Theorem 7.1 stiU holds in this case once we verify that for e > and 
for z = X or y, for any /, 

E[ka|M(|,,^l>,^V4)] — >0 asTV^oo. (4.31) 

We postpone the proof of (|4.3ip to the end of the proof. Assuming that (|4.3ip holds true, we obtain 
the CLT theorem 7.1 ( f B- Ya2) ) : the Hermitian matrix Zn — {ZN{p,q)) of size kj defined by 

converges to an Hermitian Gaussian matrix G. The Laplace transform of G (considered as a vector 
of C^, that is of {Gpq, 1 < p < q < kj}) is given for any c € by 

E[exp(c^G)] = exp[ic^Bc] 
where the K x K matrix B = {B{1, 1')) is given by: B — limjv Bi{N) + B2 + B^ with 

Bi{N){l,l') - uj{E[xnynxi,m'i]- P{1)P{1')), 
B2{lX) = {e-uo)E[xayi,i]&[xviyii] 
B3{l,l') = {t ~ uj)E[xiixi>i]E[yiiyi,i] 

and the coefficients uj,9,t are defined in Theorem 13.21 Here A = G{pg-) so that oj = 1/0| and 
6* = 1/(6*1 - ■) (ggg ^.jje Appendix). 
^From Lemma 1121 

B2{1, I') = {0 - Uj)cr^5p,qi5p' ^ (0 ~ Uj)a'^lp=q=p'=q>. 

Moreover in the complex case, B3 = and in the real case, 

BsilJ') = {e-oj)a%^i,. 

^From 4., in the real case, 

lim Bi{N){l,l') = Si,i,Loa^{l + 5p,q), 

and in the complex case, 

lim Bi{N){l,l') =Si,vLoaHp,q. 

N^oo 

It follows that 5 is a diagonal matrix given by: 

B{l,l) = (1 + 5p^q)9cr'^ = (1 + Sp^q) „f_ 2 in the real case, 
B{1, 1) = 5p,q9(j ~ 5p,q g f_^ 2 in the complex case. 

In the real case, the matrix B is exactly the covariance of the limiting Gaussian distribution G. It 
follows that G is the distribution of the GOE(fc+o. x kj^„ , a'^ / {O"^ — cr^)). 

In the complex case, from the form of B, we can conclude that the coordinates of G are independent {B 
diagonal), Gpp has variance f^/(^j — cr^) and for p ^ q, ^e{Gpq) and Qm{Gpq) are independent with 
the same variance (since B{1, 1) — ior p q). It remains to compute the variance of ^e{Gpq). Since 
the Laplace transform of 3fie(ZAr(p, g)) and '^m{Z]\[{p, q)) can be expressed as a Laplace transform 
of ZN{p,q) and ZN{p,q), we shall apply Theorem 7.1 to {Zn{p, q), Zn{p, q)) that is to the vectors 
Xi = iUip.Uiq) and yi = iUiq^Uip) in C^. We denote by B the associated "covariance" matrix of 
size 2. The variance of 3fie(Gpq) is given by ^ limAr^oo B12 (since Bn = B22 = from the previous 
computations) with 

S12 =Bi2(l)+Bi2(2)+Si2(3) 

where here Bi2(3) = 0, 

^12(1) = ujE[\Uipf\Uiqf] u;a^ and ^12(2) = {0 - uj)E[\Uipf]E[\Uiqf] = {9- uj)a\ 
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Thus, mr(5Re(Gp,)) ^ da^ 12 ^ j {2{Q'] - a"^)). We thus obtain Lemma|ll]by using that Tr(G(pe, )) = 
[N - k)tvN-k{G{pe,)) and tvN-k{G{pe^)) ^ l/Oj. 

It remains to prove (|4.3ip . The variable ajv '■— kii|'*l(|z,i|>eAri/<i) tends to in probabihty. It is 
thus enough to prove uniform integrabihty of the sequence ajv, a sufficient condition is given by 

supjv Eia^/"'] < It s^^y to t'^^t fo'^ 1 < P < %, sup^r E[|Wip|^] < oo since the Wigner 
matrix W^y has finite sixth moment and is unitary. This proves (|4.3ip and finishes the proof of 
Lemma □ 



STEP 4: We are now in position to prove that 

(eiv(A^^_,+i(M^)),...,?7v(A^^._,+,./M^))) ^ ( Ai(T4,xfc, ),..., Afc^.(T4^.xfc,)). (4.32) 

To prove (j4.32p . our strategy wiU be indirect: we start from the matrix Vk^^N and its eigenvalues 
(Ai(Vfcj^Ar); 1 < i < kj) and we will reverse the previous reasoning to raise to the normalized eigen- 
values ^Ar(Aj, ^^^(Mjv)), 1 < i < kj. This approach works in both Cases a) and b) as we now explain. 

First, for any \ < i < kj, we define Aj^-* such that 

^n{^^n) — Ai(Vfc^-,Ar), 

that is A^^ = pg^ + X^{Vk,,N)/cB,y/N. 

Since Vfc^,jv converges in distribution towards Vfc^xfcj, Ai(T4-.Ar) also converges in distribution towards 
Ai(Vfe^ xfcj )■ Hence (^iq{h^l^) converges in distribution and A^y* converges in probability towards pg-. 
Let X^j^^ = Xkj^Ni^pf) = Vkj.N +op(l) as defined in Proposition |4]3l This fit choice of A^^-* gives that 

A»(^£^) = Cjv(A^^) + e„ with e, = op(l). 

Hence, $,N{^pf) is an eigenvalue oi xj!^ — eihj- 

According to Propositions 14.11 and 14.31 on an event Ct n whose probability goes to one as N goes to 
infinity, there exists some h such that 

Aj;) = A,, (Mn - -^diag(/fe^. , ON-k, )) . 

The following lines hold on 57 jv- By using Weyl's inequalities (Lemma 13. 2p . one has for all i G 
{1, . . . ,kj} that 

XiJMn - -^diag(4^.,0^_fcj) - Ai,(Mjv) < 



N 

We then deduce that 

(^jv(Ai,(MAr)),...,CA,(A^(MA,))) = (x^{Vk^,N),---.Xu,{Vk^,N)) +op(l) (4.33) 

and thus 

(C^(A,,(M^)),...,eAr(A^(MA.))) ^ (Ai(T4^xfc,),...,Afc^(T4^.xfc,)). (4.34) 
Now, to get (|4.32p . it is sufficient to prove that 

^ li = kj-i + i; i = 1, . . . , A:j ) — > 1, as — > oo. (4.35) 



Indeed, one can notice that on the event {h — kj-i + i; i — 1, . . . , kj} the following equality holds 
true 

(eA'(\_,+i(Mw))---^A^(\_,+fc,(MA,))) = (eAr(A,,(M^)),...,C^(A^(M^))) . (4.36) 
Hence, if is satisfied then combined with imply ((i?^ . 

We turn now to the proof of (|4.35p . The key point is to notice that the kj eigenvalues of V^^xfej 
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have a joint density. This fact is weU-known if Vfe xfe is a matrix from the GU(0)E and so when 
Kj is infinite (Case b)). When Kj is finite (Case a)) and independent of N, we call on the following 
arguments. One can decompose the matrix U'^ ,^i.,HKjUKj^kj appearing in the definition (|2.7p of 
Vk^^ in the following way 

UK.xkj Hk, Uk^x kj =Qk, + 

with distributed as GU(0)E (using the fact that U'^ ,^^XJKjXkj = Ikj) and Qkj independent from 
Hkj. Hence, the law of Vk^xkj is that of the sum of two random independent matrices: the first one 
being the matrix H^j distributed as GU(0)E associated to a Gaussian measure with some variance r 
and the second one being a matrix Zk^ of the form U'^.^f.WKjUKjXkj + Qkj- Using the density of 
the GU(0)E matrix Hkj with respect to the Lebesgue measure dM on Hermitian (resp. symmetric) 
matrices, decomposing dM on Uat x (IR''^)< (denoting by \5n the unitary (resp. orthogonal) group), 
one can easily see that the distribution of the eigenvalues of Hkj + Zkj is absolutely continuous with 
respect to the Lebesgue measure dX on M" with a density given by: 

/(Ai, . . . , Ajv) = exp(-^ ^ A^) [](A, - A,)tE ("exp Tr Z^ | /((Ai , . . . , Aat), Zk, )\ dX 

i=i i<j \ L J / 

where /((Ai, . . . , Aat), ZfeJ = /exp (:^A^Tr(U'diag(Ai, . . . , AAr)?7*ZfeJ) m(rf[/) denoting by m the 
Haar measure on the unitary (resp. orthogonal) group. 

Thus, we deduce that the kj eigenvalues of VkjXkj are distinct (with probabihty one). Using Port- 
manteau's Lemma with (|4.34[) then implies that the event 

'■— \ ^n{\i (Mat)) > • • • > ^^r (A^.(M^))}f| 

is such that lim^r P(f2^) = f . By Theorem ll.il we notice that the event 

n'^ := {a^^_^ (M^) >po^+5> Xi, (M^)} f] fl [K > Pe, - ^ > A^.^_^^,^.^,(Ma,)} 

also satisfies limAr P(r2'^) = 1, for 5 small enough. This leads to (|4.35p since Cl'j^ d {U ^ i + fcj-i, i = 
1 , . . . , A; j } . 

The proof of Theorems 12.11 and 12.21 is complete. 



5 Appendix 

We recall the CLT for the empirical distribution of a Wigner matrix. 

Theorem 5.1. (Theorem 1.1 in W-Yalf ) Let f be an analytic function on an open set of the complex 
plane including [—2a, 2a]. If the entries {{WN)ii)i<i<i<N of a general Wigner matrix Wn of variance 
a^ satisfy the conditions 



(i) for I ^ I, ¥.{\{WN)a\'') ^ const, 

(n) foranyTj>0, hniAr^+oo E.*,i IE l(W^A')^;|'* l'{|(w„).,|>r,y]v} 



0, 



then N (^tr N if {^Wn))~ J fdtisc) converges in distribution towards a Gaussian variable, where figc 
is the semicircle distribution of variance a^ . 

We now prove some convergence results of the resolvent G used in the previous proofs. 
Let I < j < J+a and k such that iV — fc — > oo. 

Lemma 5.1. Each of the following convergence holds in probability as N oo; 
i) Vn (tTN-k G{pej) ~ l/0j) — ^ 0, 
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Proof of Lemma I5.H We denote by G the resolvent of the non-Deformed Wigner matrix 

WN-k/VN. 

i) By Theorem 15.11 one knows that y/N ^tr^v-fc G{p0-) — J ^^^^4^) converges in probability towards 
0. Now, we have J '^p^'^^^ = j- (see |H-P) p. 94). It is thus enough to show that 

ivN-kG{pe^)-ivN-kG{pe^) = op(l/ViV). 

Let then U^-k '■= U (resp. D^-k) be a unitary (resp. diagonal) matrix such that A^-k = U* D^^kU . 
Then, one has 

I trAT-fc G{pe^ ) - trAr_fe G{p0^ ) | = | trjv-fc {G{pg. )AN-kG{pg. )) | 

= |trjv-fc {DN-kU*G{pe^)G{pe^)U)\ 

:= \trN-k{DN^kA{pe,))\ < (r/(iV - fc))||i^jv-fe|| |1 A(pe,)|| 

where r is the finite rank of the perturbed matrix A^-k- 

One has ||DAr_fc|| < \\An\\ := c (with c = max(0i, \9j\) independent from N). Moreover on the event 

ilN := f^^^ n {\\WN-k/\^\\ < 2a + S}, \\A{pe,)\\ < [pe, -2a- (use Lemma[31]) so that we 
deduce that 

I trA,_fe(G(p,J) - trAr_fc(G(pe,))|llo„ < jf^ {pe, -2a- S)-^ ^ 0. 

Using Theorem 15. II and the fact that P(il7v) — > 1, we obtain the announced result. 

ii) It is sufficient to show that trAr_fc &{pg. ) — tr^v-fc G^{pg ) ^ in probability since, by Theorem 
15.11 one knows that tr^v-fc G^{pg) converges in probability towards f 7 — \_ y^ dpscix). 

Using the fact that Tr(BC) = Tr(CS), it is not hard to see that 

trN-kG^{p0,)-trN-kG^ipg.) = tr^-fc ((G(peJ + G(pe^ ))(G(pe^ ) - G(peJ)) 

= tr^_fc (Gipg^)Aj,^k{Gipg^) + G(pe,))G(p6*,)) 
= trAT-fc (^DN-kUGipe^)iGipg^) + G{pgyG{pe^)U*) 

:= tiN-k {DN-kA'{P6j)) 

where the matrices Diq-k and U have been defined in i). We then conclude in a similar way as before 
since on the event Oat, ||A'(pe3)|| < 2{pg. —2a — S)~'^ . 

For point iii), we refer the reader |C-D-Fj . Indeed, it was shown in Section 5.2 of |C-D-F| that the 
announced convergence holds in the case fc = 1 and for G instead of G. It is easy to adapt the 
arguments of [C-D-F| which mainly follow from the fact that, for any z g C such that 5m(z) > 0, 
Iv^ Si!=T''(^(^)")^ converges towards 5^(2). But this latter convergence was proved in Section 4.1.4 
of |C-D-Fj . □ 
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